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ANALYSIS OF AN ELECTROHYDRAULIC AIRCRAFT CONTROL -SURFACE 

SERVO AND COMPARISON WITH TEST RESULTS 

John W. Edwards 
Flight Research Center 


INTRODUCTION 


During the past decade, the use of electrohydraulic servomechanisms in aircraft 
flight-control systems has become commonplace. They are used to provide power boost, 
stability augmentation control inputs, or primary control of the vehicle. Because of 
weight and space limitations, there is always a design trade-off in the installation of 
such a system concerning the structural stiffness of its components. In a well-designed 
electrohydraulic control -surface system, a well-behaved system with a first-order- 
lag dynamic response may be achieved. The system time constant is determined by 
the position feedback gain. 

In some instances this desirable response is not achieved because excessive struc- 
tural compliance causes the system to be highly resonant. This was true of a rudder 
control servo installed in the NASA general purpose airborne simulator (GPAS). 

Attempts to damp such systems have led to the use of dynamic pressure feedback (DPF), 
in which high-passed pressure feedback provides the necessary damping. 

The literature on electrohydraulic systems is extensive (refs. 1 to 8, for example); 
however, there does not seem to be any readily available analysis of this type of sys- 
tem, which is often encountered in aircraft control -surface systems. The purpose of 
this report is to analyze a lumped, two-mass model of an aircraft control-surface 
servo and to illustrate the use of DPF. The effect of hinge -moment loading of the con- 
trol surface is included. The model is linearized, transfer functions of the pertinent 
variables are derived, and the effectiveness of DPF in damping the load resonance is 
investigated. The nonlinear and linear responses of the model are compared with 
experimental measurements of the highly resonant rudder control servo installed in 
the GPAS. 


SYMBOLS 


A 

B 

Be 

B v 


piston area, in^ 


general coefficient matrix 
coulomb friction coefficient, lb 
viscous friction coefficient, lb/ft-sec~l 



C h 

c h 6 » c h 0 

c 
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F 

f 

G,GpG 2 

H s 

I 

j = V I i 
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Kc 


K r 




K c 


K„ 


Kx 

Ki,K 2 

L 


L 


g 


l 


m p 

P F 

P L 


control- surface hinge -moment coefficient 
coefficients in hinge -moment equation 

control -surface reference chord, ft 

servovalve command voltage, V 

general coefficient matrix 

frequency, Hz 

feedback gain matrices 

control -surface hinge moment, ft— lb 

control -surface moment of inertia, slug-ft2 

spring constant, ft -lb/ rad 
servovalve flow-pressure gain, in3-sec~Vpsi 
pressure feedback gain, V/psi 
servovalve flow gain, in3-sec“Vin. 
servovalve first-stagje gain, in/V 

servovalve proportionality constant, in^-sec - VinVpsi 
position feedback gain, V/ft 

coefficients in hinge-moment equation, ft-lb/rad 
6 r 

gain constant of — (s) transfer function 

r 

servovalve coil inductance, H 

scale length of turbulence, ft 

moment arm, ft 

piston mass, slugs 

output of high -pass filter, psi 

differential load pressure, psid 
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P 1’ P 2’ P 3 

Q L 

Qi,Q 2 »Q 3 »Q 4 


R 

S 

S 

u 

V 

V a 

v g 

X 

*p 


s 


"1 

Z 1 


hydraulic source pressure, psi 
actuator chamber pressures, psi 

real poles of transfer functions 
hydraulic fluid flow rate, in 3 /sec 
hydraulic fluid flows in servovalve, in 3 /sec 
servovalve coil resistance, ohms 
control -surface reference area, ft 3 
Laplace transform variable 
control input vector 

volume of fluid between servovalve and actuator, in 3 

aircraft velocity, ft/ sec 

lateral gust velocity, ft/sec 

state vector 

piston displacement, ft 

servovalve spool displacement, in. 

auxiliary variable introduced in model 

real zero of transfer function 


P 

Pe 

A(s),A! 

6 

6 r 

^c’^t’^’^p 

9 

P 


aircraft angle of sideslip, rad or deg 

bulk modulus of elasticity of hydraulic fluid, psi 

denominator of transfer function 
control -surface deflection, rad 
rudder deflection, rad or deg 
damping ratios 

angle of incidence at control surface, rad 
atmospheric density, slugs/ft 3 
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'g 


#Vg( W ) 

(jO 

w s ,w c ,u h» ^p» w 2 

Subscripts: 

c 

g 

h 

max 

ss 


real part of complex number 

root -mean -square value of Vg, ft/ sec 

time constant of high -pass filter, sec 

time constant of first stage of servovalve, sec 

power spectrum of lateral gust 

frequency, rad/ sec 

open-loop natural frequencies, rad/sec 

command signal 
gust 

mode caused by hinge-moment loading 
maximum operating value 
steady state 


A single dot over a quantity indicates the first derivative with respect to time; 
double dots indicate the second derivative with respect to time. 


DEVELOPMENT OF MODEL 


The system analyzed in this report is an aircraft control-surface control system 
powered by an electrohydraulic servovalve (figs. 1(a) to 1(c)). The system is similar 
to the stability augmenters or power boost units in many existing systems. The mech- 
anical linkages to a pilot’s controller are not shown in figure 1(a) because only the 
electrohydraulic control system is analyzed. 

The command voltage to the servo system may be derived from the measurement 
of rate or attitude gyros or from the pilot' s controller position. The electrical input 
to the servovalve results in hydraulic fluid flow to one side of the actuator or the other 
and causes piston and control -surface motion. To make the system stable, a meas- 
urement of the piston position is fed back and summed with the command voltage. In 
some systems, this position loop closure results in lightly damped oscillations. If 
this occurs, it has been found to be advantageous to provide an additional feedback 
signal, that of the hydraulic differential pressure, across the piston. When this second 
feedback signal is properly conditioned, a well-damped system usually results. 

Both the position and pressure measurements are shown as if they were made with 
electrical transducers. This is only for ease of illustration. Systems are in use 
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which provide these feedbacks mechanically. Both methods of feedback, however, 
result in similar system responses. 

Figure 1(b) shows the mechanical components of the control-surface control system, 
a two-mass, spring-coupled system. The load moment of inertia is I and the mass of 
the actuator piston is mp. The piston is driven by hydraulic pressure, Pp,, which acts 

on the piston head area. A, and applies torque to the load through a moment arm, l. The 
linear piston position is shown as x p , and the angular surface position is indicated by 

6. There is viscous friction, By, and coulomb friction, B c , between the piston and 
the cylinder. The surface hinge moment is shown as H g . It is assumed that the com- 
pliance in the system is in the torque tube where the piston is attached, giving the 
rotary spring constant, K. For convenience, structural damping of the load is not 
considered, because structural damping ratios are characteristically small. 

Systems of the type shown in figure 1 are discussed in references 1 to 5. Ref- 
erence 1 derives pertinent transfer functions for a simplified model (m p is assumed 

to be zero and no actuator friction is included). References 2 and 3 describe a similar 
system in a spacecraft engine servo and illustrate nonlinear techniques for studying 
possible limit cycles. Reference 4 describes the aileron actuation system used on a 
B-52 airplane and mentions high -frequency cont rol -valve -gene rated noise (audible noise 
termed "honk"), an effect which will be discussed later. This report is a summary 
of many of the practical problems which may be encountered in a hydraulic system. 
Reference 5 analyzes the various configurations of aircraft hydraulic systems in detail. 


Electrohydraulic Servovalve 

A critically centered threeland, four-way, two-stage servovalve is assumed to be 
the power control device for the hydraulic system analyzed. The servovalve consists 
of a first-stage electrical torque motor, which controls pressure on the ends of the 
spool (the second stage, fig. 1(c)). The actuator used in this analysis is assumed to be 
double-ended; that is, both sides of the piston head have equal areas. Movement of 
the spool opens orifices connected to a pressure supply and sink and to both sides of 
the main actuator piston. The input to the servovalve is an electrical current, or 
voltage, and the output is the differential load pressure developed across the piston. 

An intermediate variable is the spool position, x g . Servovalve dynamics are analyzed 

in references 1, 6, 7, and 8. The succeeding development follows that in reference 6. 

To develop a model of the servovalve, the relationship of the spool position, x s , to 
the load pressure, P^, must be known. It may be obtained by apply ing the orifice flow 

equation (ref. 6, p. 41) to each of the four orifices shown in figure 1(c). The hydraulic 
fluid flows through these orifices are labeled through Q 4 . The load flow, Q L , 

satisfies the relations 
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Ql — - Q4 


_ ^3 " ^2 

The orifice flow equation shows the relationship between the fluid flow through an orifice, 
the orifice opening, and the pressure drop across the orifice. In the development of 
the orifice equation, it was assumed that the flow through the orifice was incompressible. 
Applying the orifice equation to each of the four orifices and using the relations given 
for Ql and Pl, 


Ql 




(i) 


where K v is a proportionality constant. Later in this report the effect of fluid com- 
pressibility in the actuator chamber is added to the model. This does not invalidate the 
use of the orifice flow equation here, however, because it is applied only to flows in the 
immediate vicinity of the orifice (within 1 orifice diameter of the orifice) where the 
assumption of incompressibility is valid. 

Equation (1) may be expanded in a Taylor’s series about a nominal operating con- 
dition, giving to first order the equation 


AQ l = K q Ax s - K c AP l 

where A indicates a deviation from the nominal and 


(2) 


Kq = 
K c = 


9Ql 

9x s 



= K V\/ P s-Ffj P L «*> 

K vl* S l 

= W P s-i|j P L < 3b > 


Common nominal operating conditions are Pl = 0 and x s = x s max * giving 


K q = K v> /p7 (4a) 

K v x s 

1 v b max 

K c = ^ F* (4b) 

The quantity K is the valve flow gain and is important in determining the servo posi- 
tion loop gain. The quantity is the valve flow-pressure gain. Equations (3a) and 
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(3b) are nonlinear equations for the coefficients of equation (2), whereas equations 
(4a) and (4b) provide linear approximations to the coefficients. An expression of inter- 
mediate accuracy can be obtained by expanding the square root in equation (1) as 

Ql - K v x 3 P s 2 Jl S K v x s (p s 2 - | 

If Kq and Kg are determined from this expression and the same nominal operating 
condition is used, 


<5a) 

. Kvfcsl 

K C - \ < 5b ) 


Equations (5a) and (5b) are used later in an analog simulation of the model, and 
equations (4a) and (4b) are used for linear analytical calculations. Henceforth, the A 
terms will be omitted from equation (2), giving 


QL = KqX g - K C P L 


( 6 ) 


Aircraft hydraulic systems commonly utilize a 3000-psi source pressure. The load 
pressures encountered are usually much lower than this source pressure. In this 
instance the differences of equations (3a) and (3b) , (4a) and (4b) , and (5a) and (5b) are 
small, and the use of the linear equations is justified. 

The continuity equation is used to relate the flow through the servovalve to the 
mechanical motion of the system. Mass flow through the servovalve, given by equation 
(6), is equated to flow to the actuator. The flow to the actuator is the sum of piston 
displacement and flow due to fluid compressibility (ref. 6): 

Ql = 4^ L + 12A *P (7) 

where Xp is piston displacement, V is the total volume between the servovalve and 

piston, and fi e is the bulk modulus of elasticity of the hydraulic fluid. The numeral 

12 occurs because valve coefficients are expressed in terms of cubic inches per unit. 
Equating equation (6) to equation (7), 

P L + 12Ax p = KqXg - K<,P L 

P L = ^(-KcP L - 12Ax p + KqXg) (8) 
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It is common to include the effects of leakage across the piston head and leakage 
to drain in equation (7), giving rise to leakage coefficients in equation (8). These terms 
are proportional to P L , so K c may be interpreted as a composite coefficient. Leakage 

is not considered explicitly in this report. 

In equation (8), care must be taken to insure that all units are correct. It is 
assumed that the valve coefficients are given in cubic inches per unit, P L and /3 e 

in pounds per square inch, V in cubic inches, and Xp in feet. 

The first-stage valve dynamics relating input voltage to spool displacement are 
also required for the model. The torque motor coils form a resistance-inductance 
(RL) electrical circuit which, when driven by a voltage source, exhibits a first-order- 
lag characteristic with time constant L/R and gain l/R. The torque motor is loaded 
by the first-stage hydraulic flow and the spool. Analysis of this system is complicated 
(ref. 6). Fortunately, the response of the servovalve is extremely fast and can be 
modeled adequately for low-frequency approximation as a first-order lag (ref. 8). The 
time constant, r s , is approximated from manufacturer’s test data as the 45 phase 

shift frequency of the valve frequency response. 


The response of the valve is then given by 


x s (s) 


K 


s 

1 + T S S 


e(s) 


(9) 


where e(s) is the Laplace transformed voltage input. 

The first stage of the servovalve exhibits several nonlinearities which may be 
important to the system response. The magnetic circuit of the torque motor coils will 
have some hysteresis. The magnitude of the hysteresis is given by manufacturer’ s 
test data as a percentage of maximum input current. Because the size of the hysteresis 
loop shrinks as the input signal amplitude decreases, this effect is usually not critical 
and will not be considered. Another nonlinearity is the friction between the spool and 
its housing, which gives rise to backlash in the output of the first stage of the servo- 
valve. Thus there will be a threshold value of input current below which no change in 
valve output is noted. In well-designed valves, backlash can also be ignored. 


Piston and Load Dynamics 

The equations of the piston and load dynamics are also necessary for the model. 
From figure 1(b), the equation of motion of the piston mass is 


m x 
P p 


- B v x p + F c + AP l 


(10) 


where A is the piston area. The equation of motion of the control surface is 


16 = - K 6 - 


+ H< 


(ID 
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The coulomb friction is given by- 
may be approximated as 



The aerodynamic hinge moment, H s , 


H s = |pV a 2 ScC h (ref. 9) 


where 


Ch - Ch qO + 

and 9 indicates the angle of incidence of the control surface. Thus 

H s = K x 9 - k 2 6 

where 


Kl = |pV a 2 ScC h0 
K 2 = " |pV a 2 ScC h6 


( 12 ) 


Control Law 

To complete the model the control law must be specified. It is assumed that there 
will be position feedback and high-passed pressure feedback plus a command input to 
form the voltage signal, e. The pressure feedback is modeled as 

p F (s) = TIt p s P L< S > < 13 ) 

where P F indicates a filtered pressure signal. It is assumed that the position feedback 
is derived from a measurement of the piston position, Xp. Thus the servovalve command 
voltage is 

e = K x( x p c " x p) " K p p F (14) 

where K x is the position feedback gain, Xp c is the command signal, and Kp is the 
pressure feedback gain. 

The complete model given by equations (8) to (14) is summarized in equations 
(15a) to (15g). An auxiliary state variable, Xp has been introduced to form the high- 

passed pressure feedback signal to facilitate later analysis. 
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16 = - K [5 - - k 2 6 

(15b) 

4 0e / . \ 

P L = — V~ KcP L ' 12Ax p + K q x sy 

(15c) 

1 K s 
*a = - ^ X s + 7; e 

(15d) 

X 1 = " ^ X 1 + P L 

(15e) 

P F = P L - ^ X 1 
P 

(15f) 

e = K x (x pc - x p ) - KpP F 

(15g) 


DERIVATION OF TRANSFER FUNCTIONS AND 
APPROXIMATE LITERAL FACTORS 


The model equations (eqs. (15a) to (15g)) can be manipulated to form the system 
open-loop transfer functions. Approximate literal factors of these functions are given 
in this section. To derive the open-loop transfer functions, linear equations are 
necessary. The nonlinear terms in the model are the coulomb friction, B c , and the 

valve coefficients. Linear approximations for the valve coefficients were given in 
equations (4a) and (4b) and are assumed in the following analysis. To eliminate B c 

it is assumed that a valid model may be retained by setting B c equal to zero and 
increasing By to provide the required damping at a specific frequency. It is demon- 
strated later that this is a valid technique. 

Equations (15a) to (15d) describe the open-loop system. Incorporating the above 
assumptions, these equations are Laplace transformed and written in matrix form as 


r 2 r> K 

K 


-1 


r- n 


p 



mpS^ + B y s + j 2 

l 

-A 

0 


x p 


0 

0 


K 

9 K +K 2 




6 



*1 

I 


" ll 

+ — 

0 

0 



0 

"e~ 

4 8/3 e A 

0 


■Me^q 



= 




— r s 

s+ V 

V 


P L 


0 

0 


0 

0 

0 

r s s+ 1 _ 


_*s_ 


_ K s 

0 
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The determinant of the matrix of coefficients is expanded to give the characteristic 
equation, and Cramer' s rule is used to determine the numerators of the transfer 
functions. Only transfer functions of response due to command voltage input are given 
at this time. Open-loop 6 transfer functions are of little interest because 9 inputs 
(hinge moments) represent disturbances. Closed-loop transfer functions of response 
due to 6 inputs are considered later. 


e< s >=AT5>-^V^Y + ii P) 


= JL_ 4 fte AK s K q ( 2 


(16) 


6 1 4 ^KK s K q A 

e (S ' A(s) I IV 


_ 1 'MeKgKq 

e {s> ~ A(s) V 


mps 4 + B v s 3 + 


K ( K + K 2) m p"| 9 By K 

_ + v _. 7 P I g2 + ; ( K + k 2 > s + — 


where 


(17) 
KK2_ 

J 2 

(18) 


A(s) - ( V + 1) I m p s5 + (^p> + B V ) S 4 + + Me 

K 48^ e A 2 "] [~ 4 /3 e K c m p B v ,_ _ 4/SpKK^ 

+ Z 2 + I s3 + l^v (K + K 2) + T (K + K 2) + 


KK 2 4 /3 e K c B v 


48/3 e A 2 


II 2 


+ " IV < K + K 2> + -#“ (K + K 2> s + 


4 


4^ e KKcK 2 

IV l 2 


(19) 


The above transfer functions are exact for the linear model. Thus, for a specific 
design, the parameters could be assigned values, and equations (16) to (19) could be 
used to determine the system response. This could be done easily on a digital comput- 
er. However, much insight into the dynamics of the system may be gained if literal 
expressions for various special cases can be derived. If the restating literal expres- 
sions can be shown to be valid under the assumptions of the derivations, they may be 
useful in system design and analysis. 


For the remainder of this report it is assumed that 


m P <<^2* 


In the following 


sections approximate literal transfer functions are derived from equations (16) to (19). 
The derivations begin with a basic system and become successively more complex. 

Ibe basic system is obtained by assuming that the piston mass and friction are negli- 
gible and that there is no hinge moment. This is the system analyzed in reference 1. 


Case I: Basic System 

In Case I, the basic system, it is assumed that there is no hinge-moment loading, 
that the mass of the piston is negligible compared to the load mass, and that friction 
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forces are negligible. Incorporating these assumptions in equation (19), 

i2\ 4/3 e KK c 48/3 A 2 K 


A(s) = s(r s s + 1) 
Dividing by the constant term gives 


K 48/3 e A^\ 9 

f2 + V ) S + n 2 S+ IV 


48/3 e A 2 K IV l\ 9 IK c 

A(S, - -sv— 3 <T S S + 1) + kJ S + 


S + 1 


Defining 


“h 2 = 


w; 


= 


( 20 ) 


48/3eA 

(21) 

IV 

,2 = K 
s I 

(22) 

1 l 

(23) 

= ~ 2 + ~ 2- 

co h 


K C I 

(24) 

24A 2 Z 2 


gives 


48/3 e A 2 K / s 2 2f c ^ 

A(s) = — s(r s s + l)(^ + -s + l J 


O)/ 


Hence, there is a second-order resonance in the servo response. The factor co g may 

be recognized as the natural frequency of the load coupled with the torque tube spring. 
The factor a>h arises because of the compressibility of the hydraulic fluid and is 

referred to as the hydraulic natural frequency. The two modes couple to give the load 
resonance, co c , and co c < o> s . The damping of the load resonance, £ c , is characteris- 
tically small (< 0. 1). 

The transfer functions (eqs. (16) to (18)) become 



K s K q /s 2 

12A U 2 



(25a) 
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(25b) 




Ai(s) 12A l 


gL 1 lK sKq , 

e Al(s) 12A 2 £ 2 


(25c) 


A l( s ) = s(r s s + s + lj (25d) 

Figure 2 shows this system with the position and pressure feedback loops described 

p T 

by equation (14) incorporated. The — — (s) transfer function is the ratio of equation 

X P 

(25c) and equation (25a). 

Figure 3 shows the system root locus for the position loop closure and the pressure 
loop closure. Figure 3(a) shows only the position loop closure. The pole locations, 

Pi and p£> indicate the placement of the closed-loop poles at the desired value of 

position feedback gain, K x . If w c »p. and — »p,, the -£(s) transfer function 

T s 1 e 

may be approximated by 


*P/ 0 * _ KsK q i 

e O 12A s 


Then the closed-loop — — (s) transfer function becomes 

X Pc 

K x K s Kq 

Xp 12A 

x p r (S> “ K x K s Kq 

S + 12A 

and the characteristic first-order hydraulic servo response is obtained with the time 
constant determined by the position feedback gain, K x . 

If the preceding assumptions are not warranted, the closed-loop pole locations, p^ 
and p£, may be obtained by applying the root -locus -magnitude criterion (ref. 10). If 
u c - <^s* the quadratics in equation (25a) may be neglected and 

fE/ s) = KsK q 

e v ’ 12As(t s s + 1) 
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Applying the root-locus -magnitude criterion to determine the pole locations pj^ and p 2 
gives 

K x K sKq 
12Ap(T s p + 1) 


or 




K x KgKq 

12At s 


and 

P l,2 

Thus if the quadratic terms may be neglected (i.e. , w c = w g ) the closed-loop transfer 
function following the position loop closure is 

*p KxKsKg 

x pc (S) 12 At s (s + p x ) (s + p 2 ) 

If co c ~ P^» figure 3(a) indicates that the position loop closure drives the load resonance 
mode, oOq , unstable. 

Figure 3(b) indicates the pressure loop closure (following the position loop closure) 
for two different values of r p . The position loop closure has driven the combined 

resonance unstable and established the location of poles p^ and p 2 . Also, the zeros 

at io a are exactly canceled by poles for the pressure loop closure. The usefulness of 

DPF in providing damping is largely due to this fact. The smaller value of — gives 

T P 

the best damping performance but is susceptible to static loading errors. In the limit 

— = 0 a simple proportional pressure feedback would result and static load would yield 
T P 

a static position error. As — is increased, however, the efficiency of pressure feed- 

T P 

back as a damper decreases. Thus there is an optimum trade-off between the location 
of the high-pass time constant and the static loading error. Reference 9 recommends 

placing — at -I-cj c . Figure 3(b) also shows that as the pressure feedback gain 
T P 3 
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increases, the pole at — moves toward the origin. This will cause a dipole in the 

T P 

X D 1 

closed-loop — (s) transfer function at — . The residue of this dipole may be large 
x Pc T P 

enough to cause a 5° to 10° phase shift in the system frequency response at low fre- 
quencies. Such a phase shift is sometimes attributed to unknown nonlinearities when 
it may be due to the lag-lead characteristic of the dipole. The three zeros at the origin 
cause two branches of the root locus to approach the origin from the right -half plane 
at an angle of ±60°. Thus DPF will cause instability of this model at high gains. 

Figure 3 also indicates the difficulty of attempting to achieve damping of the load 
resonance poles with conventional lag-lead compensation. The proximity of the zeros 
at a> s to the load resonance poles implies that it would be difficult, if not impossible, 

to achieve an acceptable design with conventional techniques. By canceling these zeros, 
DPF eliminates this problem. Also, the cancellation is inherent in DPF, so that exact 
knowledge of the zero location is not necessary. 

Equation (25b) illustrates the necessity of using the piston position as the feedback 
signal rather than the load position, 6. The presence of the second-order zero in the 

Xp 

numerator of — (s) provides vital phase lead in the region of gain crossover which is 
6 

c 

not provided by — (s). This point is illustrated in reference 2. 


Case II: Piston Mass and Friction Included 


In Case II the restrictions on the piston mass and friction are removed but the load- 
free case is still considered. These are the conditions which exist during ground test- 
ing of an aircraft surface servo. The analysis in the next section, in which hinge- 
moment loading is included, would apply during actual flight operation. 


Rewriting equation (19) with K 2 = 0 gives 


A(s) = s(r s s + 1) j^mpS 4 + ^ /3e y C P + B v js3 + 

^eBv 1 ^ K 48 /3e A2 V? ^eKKcmp KB V 4(3 e KK c\ 
+ V + 72 + v / s + \ IV + I + v;2 / 


+ 


Z 2 

^4/3 e KK c Bv 48/3 e A 2 K* 
+ 


IV 


IV 


VZ 2 / 
(26) 


The development of approximate literal factors of the system transfer functions for 
Cases II and HI is given in appendix A. It is shown that if K C B V «12A2, equation (26) 
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may be factored to give 


/s 2 2 £ c 

A(s) = s(— + — - 2 -s + X 



2 £ 2 

w 2 


s + l\ (r s s + 1 ) 


(27) 


The load resonance, u> c , is still defined by equation (23), and co 2 is a high-frequency 
mode resulting from the inclusion of the piston mass. 

x_ 5 

The numerators of (s) and -(s) are still given by equations (25a) and (25b). 

Pt 

In appendix A the numerator of (s) is shown to have the form 

G 

s(s + Z X ) (s 2 + 2 £pW p s + Up) 

The factors z^, a>p, and f p are derived in appendix A on the assumption that 
2 £p w p z i <<w p and tit - » f* 

l m p 


In summary, the open-loop transfer functions are: 

x P,_v 1 K s K q / s 2 \ 
T< s >- A 2(3 ) 12A ^T + 1 J 

1 K » K 1 
e l&; A2(s) 12A l 


e A 2 (s) 12 A 2 Z 2 


P L_ 1 IK s Kq / g 2 2£ 

S(s + — + — -X-s + 1 

J p 


OJZ w p 


2 £« 


2 f 2 


A 2 (s) = s(t s s + 1 ) — 5 - + — s + 1 — - + — - s + 1 




c 0, 


co 2 
<2 


co 2 


Byl‘ 

z i" T“ 


w p" ,2 


K 


“P 


(28a) 

(28b) 

(28c) 

(28d) 

(29) 

(30) 
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By 


( 31 ) 


£ 


P 


2m r 




The inclusion of the piston mass introduces an additional high-frequency mode in 
the response. This mode is the result of the interaction of the piston and the torque tube 
spring. The damping of the mode, ££> proportional to B y . 


The effect of the dynamic pressure feedback will be modified because of the inclu- 
sion of friction. The location of the zero at z 1 has a marked effect on the pressure 

loop closure, as illustrated in figure 4. The position loop has already been closed, 
and the high-frequency complex poles and zeros at cl >2 a nd Wp are neglected. The 

movement of the poles at u>z is negligible for the values of K x and Kp used to 

stabilize the load resonance poles. For instance, in the example given later, 

co c «50 rad/ sec, c o 2 ~2000 rad/ sec, and co p ~ 1500 rad/ sec. Figure 4(a) indicates 

the closure when p^ > Zp and figure 4(b) indicates the closure when p^ < z^. (Both 
Pl and z^ may be near — , so this pole location must be considered also.) 


Figure 4(a) indicates a damping action similar to that shown in figure 3(b). Figure 
4(b), however, shows the load resonance to be much better damped for the same value 
of Kp, and at larger values of pressure feedback gain, Kp, a low-frequency oscillatory 

mode emerges. It is suspected that this is the type of response mentioned in reference 
8. Furthermore, the high-frequency piston mode, » I s probably the cause of the 

high-frequency noise also observed in reference 8. There it was stated that the only 
effective means of eliminating the noise was to introduce friction on the piston. This 
observation agrees with equation (31), in which the damping of the mode is shown to be 
proportional to the friction. 

x p 

Figure 4 indicates that there will again be a dipole in the — — (s) transfer function 

x Pc 

at — , which may cause low-frequency phase lag. Also, the two zeros at the origin 
T P 

cause two branches of the root locus to approach the origin at ±90°. Thus DPF will 
not cause system instability at high gain as for Case I. 


Case HI: Hinge Moment Included 

Open-loop transfer functions. — In Case in hinge moment is included in the model, 
resulting in equation (19), the characteristic equation. In appendix A it is shown that if 

40 e K c B v KKo 48/? e A 2 

IV (K + K 2) + T7F « -^^< K + K 2> 
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and pg « oj c , equation (19) may be factored to give 


A(s) = (t s s + 1) (s +p 3 ) 


, w ch 


2 ^ h A 

+ s + 1 

w ch / 


K w 2h" 


2S2h J 

+ s + 1 

w 2h 


where the subscript h indicates the effect of hinge moment. It is noted in appendix A 
that 


w ch ~ 
^ch = 


/K + K 2 
iT" 

J k + k 2 


w c 




hence the stability and damping of the load resonance poles of the loaded system will be 
slightly improved. It is also noted that o>2h = w 2 a *id £ 2 h = ^2> indicating the inde- 
pendence of the piston mode from hinge-moment effects. 

x p 6 P L 

The numerators of -^(s) and -(s) remain unchanged, but (s) must be 

G 6 0 

p T 

reconsidered. Appendix A indicates that the numerator of — - (s) may be approximated 

0 

by 


( s2 + 2 £ph w phs + l)(s 2 + 2£ t cu t s + cut) 

k + K 2 K 

if “I— - « — 2 — and cut « w ph* 

l nip 


In the preceding expression 


w ph — w p 


£ph ~ £p 



, 1 r± 

£ t 2 I \ K 2 
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The factored transfer functions for Case m are given in table 1. The simplified 
models of Cases I and H (eqs. (25a) to (25d) and (28a) to (28d)) may be obtained by 
setting mp, B v , or K 2 equal to zero. The subscript h has been eliminated from the 

parameters in table 1. 

Figure 5 shows the pressure loop closure (following the position loop closure) for 
Case in. The load resonance poles at a> c will close on the zeros at cot* Thus the 

location of these zeros is critical to the stability of the system. If they lie in the region 
of co c , DPF will be ineffective as a damper. Table 1 shows that is proportional 

to the square root of the hinge moment, K 2 , and that £ t is proportional to By. Thus 

some amount of friction on the piston is desirable, because it causes the zeros at 

to move toward a region of higher damping. Of course, piston friction is also needed 
to supply damping for the high-frequency mode at w 2 * Also, the low-frequency oscil- 

latoiy mode which occurs at high DPF gains and low values of position loop gains 
(fig. 4(b)) does not occur with hinge-moment loading. 

Hinge-m oment loading errors. -Only open-loop transfer functions of system response 
due to a command voltage input have been considered. These transfer functions and 
their approximate literal factors, which have been developed, should make it possible 
to design a system with a suitable closed-loop response. Once this has been accom- 
plished, the sensitivity of the closed-loop system to hinge-moment loading is of interest. 

To stucfy this effect, the — (s) transfer function of the closed-loop system will be 

determined. Thus the open-loop matrix used to derive the open-loop transfer functions 
must be augmented by the control law (eqs. (15e) to (15g)). The voltage signal is now 


e - K x (x Pc - x p ) - KpPp 
= Kx ( X P c ■ X p) “ K P P L + ~ k J) k 1 


and the x s equation becomes 


T s Xg x g + K s 


K x x P c K x x p K p P L + KpX^j 


where 


x l = -— x 1 + P L 
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Using the last two equations, the closed-loop system matrix is 


2 D K 
+ B v s + — 


K_ 

l 


-A 


K_ 

U 


48 li e A 


K + K 2 
I 




Y 


Pc 


0 0 


(32) 


Q S + 1 


'P 


-1 


0 0 


The numerator of the — (s) transfer function is 

ft 


\ (m n s 2 + B v s + 


_K 

l 2 


[(. ♦ + 1 )(. + 


(33) 


/ , 1 \ 

~48p e A 

(t s s + l) 

4£eKqK s K x 

r^) 

v s| 

1 + V 


Unfortunately, there are no simple literal factors of equation (33). It is a complex 
function of the feedback quantities K x , Kp, and Tp. It is possible to calculate the 

static closed-loop load position, 6, due to hinge-moment loading by computing only the 

constant term of the -x(s) transfer function. Equation (33) provides the numerator 

ft 

constant, and the system matrix furnishes the denominator constant: 



KKpK c + ^ 2 K 1 K q K s K x A 
KK 2 K c + 2 2 (K + K 2 ) KqK s K x A 


This expression gives the exact steady-state hinge-moment loading error. If the terms 
involving K c are neglected, 



Ki 

K + K 2 


This result could have been anticipated, because the high-pass filter in the pressure 
feedback loop isolates the load from closed-loop dynamics in the steady state. For a 
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specific problem, the system constants may be inserted into equation (33), and the re- 
sulting polynomial may be factored to give the frequency dependence of the hinge- 
moment loading error. 


COMPARISON OF GPAS RUDDER AND MODEL RESPONSE 


GPAS Rudder Electrohydraulic Servo 

Tfie GPAS is a variable -stability JetStar airplane operated by the NASA Flight 
Research Center. The basic control surfaces were disconnected from the left-hand 
pilot’ s controllers, and irreversible electrohydraulic servos were installed. Thus 
the control -surface positions may be electrically controlled for airborne simulation 
purposes. Details of the GPAS system are given in reference 11. The GPAS rudder 
is used in the following discussions to illustrate the analysis presented in the previous 
sections. Thus the surface deflection, 6, of the previous sections becomes the rudder 
deflection, 6 r , and the angle 9 is interpreted as the aircraft sideslip angle, /?. 

Figure 6 shows the location of the rudder hydraulic actuator and its relationship to 
the rudder surface. The electrohydraulic system consists of a servovalve, which con- 
trols fluid flow to the double-ended hydraulic actuator. The system operating pressure 
is 3000 psi. The rudder surface is mass -balanced. The constants required in equations 
(15a) to (15g) for the rudder model are given in table 2. Appendix B discusses the 
problems encountered in determining these constants. In order to install the servo, the 
rudder torque tube had to be cut and the moment arm linkage had to be inserted. In 
the basic JetStar there was no power boost to the rudder, and in the GPAS modification 
the rudder is permanently connected by cables to the right-hand pilot’ s rudder pedals. 

It was not necessary to include these rudder pedals (an extra mass-spring system) in 
the analysis. 

An electrical pressure transducer is used to provide the pressure feedback signal. 
Figure 6 shows that several feet of hydraulic tubing are required to connect the servo 
valve, the pressure transducer, and the actuator. In recent years it has become 
common to use an integral servo valve/actuator system with built-in mechanical position 
and pressure feedbacks. Such units are useful for large installations and for specific 
applications, because major mechanical modifications are necessary to alter the sys- 
tem characteristics. For small installations, the versatility of electrical feedback 
signals may be important. Either method will provide the same system characteristics. 
The rudder position transducer is a rotary potentiometer located beneath the actuator 
and strapped around the torque tube (fig. 7). 


Early testing of the rudder servo revealed a stair-stepping response to a step input 
command. It was recognized that the response was caused by the static friction on the 
piston. The problem was solved by superimposing a constant -amplitude, 60 -hertz 
dither signal on the command signal. The amplitude of the dither signal was adjusted 
to effectively "break" the static friction while leaving the overall response unchanged. 
Hence, only coulomb friction, B c , is included in this analysis. 
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Rudder and Model Response 


To compare the rudder response with the model response, the model equations of 
motion (eqs. (15a) to (15g)) with the constants given in table 2 were mechanized on an 
analog computer, and step and frequency responses were obtained. Two nonlinearities 
were programed: the valve flow-pressure coefficient, K c , and the coulomb friction, 

B c . Equations (4a) and (4b) or (5a) and (5b) were used for the valve coefficients, K c 

and K n . Thus a function switch selected either equation (4b) or equation (5b) as the 

term included in the pressure equation. The coulomb friction, B c , could be eliminated 

by setting a potentiometer to zero. As indicated in the section on DERIVATION OF 
TRANSFER FUNCTIONS AND APPROXIMATE LITERAL FACTORS, the viscous friction 
coefficient, By, may then be increased to yield the desired damping at a prescribed fre- 
quency. For the linear model, B y was increased from 500 lb/ ft-sec -1 to 

2000 lb/ft-sec -1 . A linear model is desirable because root loci and frequency response 
analysis techniques may then be used. Thus it is important to verify the suitability of 
the linear model. A digital program, CONTROL, which is described in appendix C, was 
used to calculate the root loci and frequency responses of the linear model. 

The step and frequency responses of the rudder and the nonlinear model are given 
in figures 8 to 10. The rudder tests were made on the ground, so there was no hinge- 
moment loading. The test input was always a 1° command. All responses were obtained 
with K x = 1090 V/ft. The pressure feedback gain, K p , and the time constant of the 

high-pass circuit, T p , were varied. Figure 8 shows the step response of the rudder 

and the model for — = 8.3, 16.6, and 33.2 rad/sec. Values of K p range from 
T P 

0.0048 V/psi to 0.0131 V/psi. At values of K p less than approximately 0.004 V/psi 

the response becomes unstable. The 60 -hertz oscillation in the rudder response in 
figure 8(a) may be caused by the dither signal. The agreement between the rudder and 
nonlinear model step responses shown in figures 8(a) to 8(c) is considered good. 

Frequency responses of the rudder, the nonlinear model, and the linear model for 
two different values of K p and r p are shown in figures 9 and 10. In figures 9(a) and 

9(b) the load resonance at approximately 7.5 hertz is evident. Note also the charac- 
teristic first-order lag low-frequency attenuation, and the notch characteristic at 
f « 12 hertz. The first-order lag characteristic can be explained as follows: The peak 
at 7.5 hertz and the notch characteristic at 12 hertz are caused by a pair of complex 
poles and zeros with very low damping ratios. These poles and zeros cannot yield the 
lower frequency phase attenuation shown in the figures. This phase shift is caused by 
the real pole, p^ (fig. 3(a)). 

In figures 10(a) and 10(b) the load resonance was damped by the larger pressure 
feedback, K p . The agreement between the rudder and model frequency responses in 

figures 9 and 10 reflects the agreement in the step responses. 
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Figures 9 and 10 also indicate that only slight differences exist between the linear 
and nonlinear models. The main difference is in the low-frequency phase shift, where 
the nonlinear model has more phase shift than the linear model. The rudder phase 
shift curve is considerably below those of both models at low frequencies. This is 
especially true with the high damping shown in figures 10(a) and 10(b). This may be 
caused by nonlinear effects not included in the model. Possible sources of the phase 
shift are torque motor hysteresis, first-stage servovalve backlash, or actuator static 
coulomb friction. 


Figures 8 to 10 indicate that the rudder servosystem is approximated well by the 
nonlinear model and that the linear model is a valid approximation of the nonlinear 
model. Thus a study of the utility of DPF with the linear model would allow more 
general conclusions to be drawn with regard to the damping action of DPF. 

The effect of hinge-moment loading of the rudder surface is shown as a final check 
of the model in figure 11. Figure 11(a) shows the rudder ground test step response, 
and figure 11(b) shows a step response recorded during flight. The loading causes the 
servo to be better damped, and the resonant frequency (in the step response) increases. 
Figures 12(a) and 12(b) show the step response of the nonlinear model with and without 
hinge-moment loading. Inspection of these responses leads to the same conclusions 
as were drawn from figure 11. 


Table 3 compares the approximate literal factors developed in the analysis section 
for Cases I, n, and III to the exact values obtained with the CONTROL program. In 
Case I the literal factors are exact. In Cases II and III the approximate factors are 
close enough to the exact values to make a first analysis possible. For the rudder 


model 


l 2 


72 slugs and m p = 0. 166 slug. 


Thus the condition 


m P <<: is satisfied. 


Table 3 also shows that the high-frequency poles at 0 J 2 are overdamped. Thus 

there should be no high-frequency control -valve -generated noise ("honk”) in the system. 
The damping is attributed to the large value of B v . Of course, B v was adjusted to 

match the system' s response at the load resonance frequency (approximately 50 rad/sec). 
Thus a separate calculation with a lower value of B v (corresponding to the higher 

frequency oj 2 ) ma Y be in order to estimate the damping of this high-frequency mode 

more accurately. 


ANALYSIS OF DYNAMIC PRESSURE FEEDBACK 


The effect of dynamic pressure feedback on system stability was analyzed, using 
the GPAS rudder servo for illustration. In addition, the effects of K x , K p , and 

on cases I, II, and III were studied. 

T P 


The CONTROL P rogram was used to obtain the root loci of the linear rudder model 
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for the three cases with K x , Kp, and as parameters. The loci are plotted in 

figure 13. Loci of the load resonance poles are plotted for three values of — and 

T P 

four values of K x . Only the third quadrant of the s plane is shown. Dots along the 
loci indicate gain increments of AKp = 0. 00348 V/psi. 

Figure 13(a) gives the loci for Case I. The GPAS system design called for a closed- 
loop rudder servo bandwidth of 50 rad/ sec. This goal is most closely achieved by the 
position feedback gain, K x , of 1090 V/ft and is the value used in the system. This 

value of K x places p^ at 42 rad/sec; however, this system would be objectionable 

because of the limited damping action of the DPF gain, K p . Figure 13(a) also illustrates 

a point made in connection with figure 3(b)— that high K p gains may drive the system 

unstable. Figure 13(a) shows that better damping performance can be achieved by re- 
ducing the position feedback gain, K x , with the attendant loss in system bandwidth. This 

model of the system (Case I) is unrealistic because friction on the piston was neglected. 

Figure 13(b) gives the loci for Case 13. The damping for K x = 1090 V/ft is reason- 
able, and, again, as K x decreases, the damping performance of the DPF improves. 
When K x = 436 V/ft and 654 V/ft, the effect is the same as that noted in figure 4(b), 



Kp increases. If Kp were increased from zero, the load resonance oscillations 

would become well damped (they would, in fact, disappear), leaving a desirable step 
response. If K p were increased still further, the low-frequency mode would appear 

in the response. 

Figures 13(a) and 13(b) indicate that for many designs there will be a value of Kp 

beyond which DPF does not increase damping; in fact, damping may decrease. More 
importantly, DPF will then also reduce the system bandwidth (through the frequency 
of the load resonance poles). 

Figure 13(c) gives the loci for Case III. The load resonance poles at w c close on 
the zeros at cjf Differences between figures 13(b) and 13(c) illustrate the difference 

that may be expected between ground tests and flight operation of aircraft hydraulic 
servos. The frequency of the load resonance poles is lower in figure 13(b) (56 rad/ sec 
to 58 rad/ sec) than in figure 13(c) (62 rad/ sec to 64 rad/ sec). This could have been 
anticipated by noting from table 1 that a> c increases as K 2 increases. In figure 13(c) 

there is a maximum damping that is practical. The loci for the different values of K x 

are similar, and all yield acceptable damping. In addition, the change in the loci of 
figure 13(c) for changes in r p is small relative to the larger changes of figure 13(b). 
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Also, increments of Kp in figure 13(c) cause smaller load resonance pole motion than 
in figure 13(b). 


The root loci of figures 13(a) to 13(c) give insight into the behavior of the load reso- 
nance poles. To augment this information, the closed -loop frequency response of 
the system indicates the system bandwidth, phase margin, and low-frequency phase 
shifts. Figures 14 and 15 give the frequency response of the linear rudder model for 


variations in K x and K p with 


held constant at — = 16.6 rad/ sec. 


The hinge- 


moment coefficient, K£, was set at zero, which corresponds to the ground test value. 

In figures 14(a) and 14(b), in which K p = 0.00348 V/psi, K x has strong control over 
the resonant peak at 7 hertz to 8 hertz. However, lowering K x also lowers the band- 
width (fig. 14(b)). In figures 15(a) and 15(b), in which K x = 1090 V/ft and K p is varied, 
K p does not exercise strong control over the magnitude of the resonance, but it does 
draw the resonance to lower frequencies (fig. 13(b)). 


Figures 16(a) and 16(b) illustrate the change in closed-loop response when hinge- 
moment loading is included. When hinge moment is included, the frequency of the 
resonance peak increases, its magnitude decreases, and the phase margin increases. 
These effects bear out the comments made in connection with figures 11 and 12. 


The following conclusions about the damping action of DPF may be drawn from 
figures 13 to 16. If w c is in the frequency region of the desired closed-loop bandwidth, 

the position loop closure will drive the load resonance poles toward the unstable region. 
With — set at approximately ^ co c , DPF will always have some damping effect and 

T p 3 

will usually be quite effective. However, if K x is large, the resulting DPF damping 
action may not be satisfactory. Decreasing K x always allows better damping action 
by DPF, but this causes a loss of bandwidth. There will usually be a value of Kp 

beyond which DPF does not increase damping (it may, in fact, decrease damping) and 
directly lowers the closed-loop bandwidth. Inclusion of hinge-moment loading decreases 
the sensitivity of the load resonance poles to variations in K x , Kp, and r p . In addi- 
tion, it raises the frequency of the load resonance, lowers the magnitude of the reso- 
nance peak, and increases the phase margin. Thus, the operation of the system will 
be more stable during flight than during ground operation. 


The use of DPF in an aircraft control-surface servosystem makes the system sensi- 
tive to variations in hinge-moment loading. Hinge-moment loads arise out of the vari- 
ations in the angle of incidence at the surface caused by aircraft maneuvers, control- 
surface deflection, or atmospheric turbulence. Whether any of these effects unduly 
compromise the servo response is considered in the following discussion. 

For the aircraft rudder servo being analyzed, load variations are attributed to 


25 



changes in sideslip angle, (3, rudder deflection, 6 r , and turbulence, y3g. The response 

of an aircraft to a rudder step is caused primarily by the excitation of the Dutch roll 
mode. Ignoring the effect of the roll and spiral modes, the response of the GPAS to a 
rudder step for a selected flight condition can be expressed as 



1.76 

s 2 + 0.265S+ 2.4 


This function was mechanized in the analog simulation, and the nonlinear rudder model 
was used to obtain the time history shown in figure 17. The frequency of the load vari- 
ation was small enough so that the rudder servo was able to maintain constant com- 
manded Xp following the initial transient; however, the rudder surface position, 6 r , 

varied as it was forced by (3 . Because the frequency of the Dutch roll oscillation was 
much less than the bandwidth of the servo, the variation in 6 r may be approximated by 

AS r “ (?) 3 A £ 


Load variations will also be caused by atmospheric turbulence. Lateral gust veloc- 

v e 

ity, Vg, induces an effective incremental sideslip of % = y and will be sensed by the 

3 . 

servosystem as a changing hinge moment. The response of the rudder to turbulence 
was simulated using the gust power spectra: 


§v g M = 



where V a is the aircraft velocity, L g is the turbulence scale factor, and o- g is the 

root -mean- square gust velocity. This power spectrum may be realized by passing 
white noise through an appropriate shaping filter. 

Figure 18 shows the GPAS rudder servo response to turbulence at a cruise flight 

L_ 

condition. A moderate turbulence level of cr„ = 5 ft/sec was used. The ratio W as 

& v a 

set at 10 seconds. The response of 6 r indicates that the response of the servo to load 
variations caused by atmospheric turbulence is acceptable. 
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CONCLUSIONS 


A lumped parameter, two-mass system driven by an electrohydraulic servo valve 
was analyzed. The results of the analysis of a model of the system were compared 
with test results from a highly resonant rudder electrohydraulic control system. The 
comparison led to the following conclusions: 

1. The model equations of motion accurately represented the test system. The 
calculated model root loci, frequency responses, and step responses agreed well with 
test results. A linear model agreed well with a nonlinear model that included the non- 
linear flow equation of the servo valve and coulomb friction on the actuator. 

2. Approximate literal transfer functions yielded accuracy acceptable for initial 
analysis if the piston mass was much less than the effective load mass. 

3. The system response was characterized by a load resonance and a piston-hydrau- 
lic mode. If there was insufficient friction on the piston, the high-frequency piston- 
hydraulic mode may have been underdamped and could have caused high-frequency 
control -valve-generated noise referred to as "honk." 

4. When the position loop was closed, the load resonance poles were driven toward 
the unstable region. By locating the high-pass time constant of the pressure feedback 
at approximately one-third of the load resonance frequency, dynamic pressure feedback 
(DPF) always provided damping for the load resonance pole. However, if the position 
feedback gain were too high, damping provided by DPF may have been insufficient. 

5. A maximum practical DPF gain usually existed. Increasing the DPF gain 
beyond this point directly lowered the system bandwidth and possibly decreased the 
damping. 

6. Operation with hinge-moment loads (corresponding to flight operations) resulted 
in a higher load resonance frequency, smaller frequency response peaking, and in- 
creased phase margin than when hinge-moment loading was not included. 

7. Hinge-moment loading reduced the sensitivity of the load resonance poles to 
variations in position and pressure feedback gains and in the high -pass time constant. 

8. The response of the rudder servo to load variations caused by atmospheric 
turbulence was acceptable. 


Flight Research Center, 

National Aeronautics and Space Administration, 
Edwards, Calif., March 3, 1972. 


27 


APPENDIX A 


DETERMINATION OF MODEL CONSTANTS 


In the determination of the constants in table 2 some difficulties were encountered. 
Although the load moment of inertia, I, and the piston mass, m p , should be relatively 

accurately determined quantities, a problem may arise in deciding the location of 
the spring, K. This location will obviously affect m p and I. If a distributed spring 

is encountered, the lumped-parameter model of this report may be completely 
inappropriate. 

The spring constant, K, should be accurately determined, because it is directly 

related to w s , the load resonance frequency. Equation (25a) indicates that — (s) has 

a notch characteristic at co s . Thus K may be determined by calculation, measure- 
ment, or (knowing I) determination of o> s . 

The friction, B v and B c , may be experimentally determined with low-frequency 
triangular wave inputs to the closed-loop system. 

It was not necessary to include the rate limit in the servo analysis. Electronic 
safety trip circuits cut off the GPAS hydraulic system before this rate limit is encoun- 
tered, and rate limiting is rarely encountered in the operation of the system. This 
reflects the fact that Pp, «Ps* 

The valve quantities K s , Kq, Kc, and r s can be obtained from manufacturer’ s 
data. The first-stage servo valve gain, K s , is a temperature -dependent quantity be- 
cause of its dependence on the coil resistance, R. The resistance of copper wire 
increases by 22 percent in rising from 68° F to 167° F. Thus, the valve gain, K s , 

which is inversely proportional to R, will decrease by 18 percent during this tempera- 
ture rise. Because K s is a factor in the position loop gain, the servo response will 

also deteriorate. This problem may be sidestepped by using current to drive the servo- 
valve. 

The bulk modulus of elasticity of the hydraulic fluid is difficult to estimate. Hand- 
book values of (3 e = 250,000 psi for typical petroleum-based fluids are usually too 

high. In reference 6 the strong dependence of /3 e on entrapped air and container 

compliance is illustrated. At the high pressures used in aircraft hydraulic systems, 
entrapped air will be dissolved soon after the system begins to operate. However, 
reference 7 indicates that entrapped air may still cause problems: Dissolved air may 
come out of solution when the fluid passes through a valve orifice due to the pressure 
drop. Reference 6 states that a (i e of 100, 000 psi yielded good results for many 

studies. The value used in this report, 0 e = 50,000 psi, was determined by attempting 
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APPENDIX A - CONCLUDED 


to match the analog simulation response to the measured response. It is believed that 
values of /3 e between these two values will give good results in analyses of hydraulic 

systems. The value used in an analysis may be chosen according to the relative amount 
of hydraulic tubing and connections contained in the system. The rudder control system 
analyzed in this report contained a relatively large amount of tubing and connections, 
and the smaller value of p e produced good results. On a close-coupled installation 

(for instance, with the servovalve mounted on the actuator) the larger value recom- 
mended in reference 6 would be appropriate. 


29 



APPENDIX B 


DEVELOPMENT OF APPROXIMATE LITERAL TRANSFER FUNCTIONS 


The approximate literal transfer functions for Cases n and III are developed in this 
appendix. 


CASE II: PISTON MASS AND FRICTION INCLUDED 


In Case II the restrictions on the piston mass and friction are removed, and the 
hinge-moment loading is neglected. 


Rewriting equation (19) with K 2 equal to zero gives 


A(s) = s(t s s + 1) 


mpS 4 + ( 4 ie^ + Bv) s 3 + (^£ + 


4/? e Bv K c K 48 feA^ 2 


V 


+ 9 + 

l 2 


v 


/4/JeKKcmp KBy 4/J e KK c \ / 4/? e KK c B v 48feA 2 K 
IV + I + vi 2 / + V iv + IV 


(36) 


In the constant term 


40 e KKcBv 48/3 e A z K 

« 


IV 


IV 


or 


K c B y « 12 A 


will usually prevail. Hence, dropping the first term of the constant and dividing through 
by the remaining term, 


48/3 e A 2 K 

A(s) = “Tv — S ( T s s + X ) 


IVm 


48/3 e A^K S + \l2A 2 K ‘ r 48/3 e A 2 K, 


/lm p Kc IVBy \ 3 


Vm p IB v K c IV n 2 + 

+ \48/3 e A 2 + 12A 2 K + 48/3 e A 2 ^2 + S 


m 


pKc ByV 


IK, 


+ 


12 A 2 48/3 e A 2 12A 2 l 2 


s + 1 


(37) 
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This expression has the form 





2f 2 > 

+ — s + 1 
w 2 j 


-c 2 w 2 


(-£- 

\o^co2 


+ 


fl£-\ s 3 

oj c w|/ 


+ 


1 1 4 “ 2^c\ 2 / 2 £c 2 -2 \ n 

—5- + + s z + l— + Is + 1 

cog 0^2 to c u)2 J \ w c w 2 / 


( 38 ) 


The last two terms of the coefficient of s 3 in equation (37) may be recognized as 


the quantity already defined as in equation (23). Then the coefficient of s 4 in 

w c 

equation (37) is 


2 ? 

w c ' m 2 = 


48/3 e A^K 

IVm n 


which may be solved for co£ (using equation (23) for co c ) as 


2,2 


2 _ 
w 2 - 


VK + 48/2 e A £ 
Vl 2 m T 


To determine £ 2 » K c « 1* (See table 2.) Hence, setting K c equal to zero in 
equation (37) and observing the coefficients of s 3 and s (where f 2 occurs) results 
in the expressions 

2 £c 2£ 2 IVBy 

2 o “ o — (39) 

w c t ‘ ,, 2 ^ 2^0 48/3 e A 2 K 

and 

2 C c 2 S2 B y V 

“C + -2 = (40) 

(Although t 2 also occurs in the coefficient of s 2 in equation (38), that coefficient is 
dominated by the other terms in the coefficient.) 
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The load resonance, w c , remains in the vicinity of the jo> axis near the zero at co s 
as shown in figure 3(a). Thus £ c «l, whereas £ 2 * Setting £ c equal to zero and 
solving equations (39) and (40) gives 


h ~ 


VBy 

96£ e A 2 



(41) 


^2 = 


VBy 

96/3 e A 2 


u> 2 


,2* K 

c I 

used to compute 


Because « - 7 - these expressions are nearly equivalent and equation (41) will be 
c I 


x p 6 

The numerators of — (s) and — (s) are still given by equations (25a) and (25b) and 

e e 

P L 

do not require further analysis. The numerator of (s) does require further analysis. 

e 


Setting Ko equal to zero in equation (18), the numerator of — (s) contains the expres- 
^ e 

sion 


O B v 2 
S ' 3 + — s* + 
m p 


K 


l 2 m r 


K\ B v K 

+ ^' 8+ tor 




(42) 


This expression has the form 


;(s + z 1 )(s 2 + 2^ p cjpS + oj p 2 ) = s|s3 + (Zi + 2f p u)p)s i: 


+ ( 2 £p w p z l + w p) s + z l w p] 


(43) 


Pt 

In Case I the ~~( s ) numerator contained s 2 . Assuming that the low-frequency 
response will not be grossly affected by setting mp ^ 0 and By ^ 0, it is anticipated 
that z i will be small and o>p will be large. Therefore, 

2 £pV 1<<C0 p 
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Comparing the coefficients of the first power of s in equations (42) and (43) results in 
the expression 



K K _ K 

T 2 — + T -72“ 

l m p i nip 


(44) 


since 


K 

;2 

l mp 



Then the constant terms of equations (42) and (43) give 



BvK 

Im p 


or 


Zi = 


B v V 


To find the coefficients of are equated: 


Z 1 + 2 £p w p = 


By 


m r 


By 


’P 2oj p \m p 1 
In summary, the open-loop transfer functions are 


\ _ ^ / l l±\ 

) 2 ^p \ m p I / 


1 K s K q 


Xp 

T (S)_ A 2 (s) 12 A \ u 


s 2 , 

~~T + 1 




1 K C K, 


s q 


e (S) A 2 (s) 12 AI 


P L \ 1 

TT( s ) = 


IK s K q 

A 2( s > 12 A 2 ^ 2 


s(s + z i) —r + 




2£, 


0) r 


-s + lj 


A 2 (s) = s(r s s + l)^ + -s + ll(- 


2£n 


2^ 




S + ^ 


(45) 


(46) 


(47 a) 


(47b) 


(47c) 


(47d) 
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where 


V 2 


2 

VK + 48/3 e A 2 Z 2 

I 


^2 = 

VZ 2 m p 

K 


, , 2 _ 

K 

£ 2 m p 


w s ' 

I 

B v / i 

iA 

o 

48/3 e A 2 1 2 

2aj p \ m p 

- T ) 

*5 = 

IV 

VB v / 

I\ 2 

1 

1 1 

48/3 e A 2 K ' 

K / w c“2 

-C 2 = 

co| + 


CASE IE: HINGE MOMENT INCLUDED 


In Case IE the restriction on the hinge moment is removed, giving equation (19), a 
characteristic equation. In the coefficient of s it is noted that usually 


^e^cBy 


IV 


(K + K 2 ) 


KK 2 48/3 e A 2 


(K + k 2 ) 


Thus, neglecting the first two terms of the coefficient of s in equation (19) and dividing 
through by the remaining term gives 


48/3 e A 2 (K + K 2 ) 


IV 


( T s s + !) 


IVm P .5 + 

IKgm p IVBv 

48/3 e A 2 (K + K 2 ) 

12A 2 (k + K 2 ) 48/lgA 2 (K +K 2 ) 


Vm r 


IB V K C 


IVK 


+ [j48/3 e A 2 + 12 A 2 (k + K 2 ) + 48/3 e A 2 z2(K + K 2 ) + K 


i-} 
+ k 2 J 


n , 

+ |_ 12 A 2 + 48/3 e A 2 + 12A 2 Z 2 (K + K 2 ) 
This expression has the form 


K C m £ 


VB v IKK C 


] 


+ s + 


kk 2 k c 


12A 2 Z 2 (K + K 2 ) 


(48) 
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+ 


2 £ch 

w ch 


s + 



+ 


2 j2h 

u) 2h 


•s + 1 



where the subscript h indicates the effect of hinge moment. 
It is obvious that 

_ KK 2 K c 
P3 12A 2 Z 2 (k+K 2 ) 


(49) 


The factor p 2 is small compared to the other factors in the expression. Thus fac- 
toring the root (s + P3) from equation (48) leaves an expression very similar to 
equation (37): 


480 g A 2 (k + K 2 ) 


IV 


(t s s + 1)(S +p 3 ) 


L Vm p c 4 , 

IK c m p 

IVBy 

48/3 e A 2 (K + K 2 ) 

12A 2 (K + K 2 ) 

+ 48/3 e A 2 (K + K 2 )_ 


Vmp IB V K C IVK 

L. — + 

48/3 e A2 12A 2 (K + K 2 ) 480 e A 2 Z 2 (K + K 2 ) 


_!_1, 
K + K 2 J 


m r 


VBv 


IKKc 


12 A 2 + 48/3 e A 2 + 12A 2 Z 2 (K + K 2 ) 


i) + I 


Because of the similarity between this expression and the expression in equation (37), 
the procedure is exactly the same as that in Case n, that is, identify 

2 K + K 2 
^sh ~~ j 

2 48/3 6 A 2 Z 2 (k + 

w hh = IVK 



1 1 1 

2 “ 9 + 9 

w ch w sn ^hr 

IK c K 

£ ch “ o 9 / \ 

24A 2 Z 2 (k +K 2 ) 


w 2h 2 


VK + 48/? e A 2 Z 2 
IVZ 2 m p 
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Note that 


£2h 


/ I \ 2 

96/? e A 2 \ K + k 2/ Wch W2h 


k’ch _ 
^ch “ 


J 


K + K 2 

r 






k + k 2 


K 


Thus the stability and damping of the load resonance will be slightly improved. Also, 
w 2h ^ 2h are same as i n Case n. The independence of and J 2 h °f the 

hinge moment reflects the fact that this is a piston mode isolated by the load mass from 
hinge-moment effects. 


Xp 5 

The numerators of -^(s) and — (s) are still given by equations (25a) and (25b), but 
P L 6 6 P L 

the numerator of (s) has changed. The (s) numerator (eq. (18)) now contains 

e g 

the factor 


4 . ?v _•} . / K 


s'* + -r~ s 3 + 

m p 


^ K + K 2 \ 9 ^ MK + K2) . ^ JEL 

l 2 m p 1 / V m p 


and has the form 


(s 2 + 2£phWphS + w p ^)(s2 + 2^^t s + ^t 2 ) 


It is anticipated that 
K + K 2 
neglecting — ; 


w ph * w p and w ph » wf 
... K 

in comparison with — 

l mp 


Equating coefficients of s 3 and again 
results in the expression 


Then 



K 

Z 2 m p 




KK 2 

U 2 m p 
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k 2 

I 


Equating coefficients of s. 


2 £ph w ph + 2 £t w t- p 


Using equation (46) for fpj^ 

„ . i nr r^v , jw.j_ 

t 2>/ k 2 [ m P ph2 " P \ m P 1 )\ 



The factored transfer functions for this case are summarized in table 1. The 
simplified models of equations (31) and (25) may be obtained by setting nip, By, or 

K 2 equal to zero. Thus the subscript h was eliminated from the parameters in table 1. 
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CONTROL: A DIGITAL PROGRAM FOR THE DETERMINATION OF 
TRANSFER FUNCTIONS, ROOT LOCI, 

AND FREQUENCY RESPONSES 


The CONTROL digital program was used to generate the transfer functions, root 
loci, and frequency responses of the linear model. The program uses the state vana 
ble formulation of the linear equations of motion: 


x = Fx + Bu 


(34) 


where 

x is the n X 1 vector 

F is the n X n matrix 

B is the n X m matrix 

u is the m X 1 vector 

n is the dimension of the system, 

m is the number of inputs into the system. 

From the closed-loop system matrix of equation (32), n = 7 and m = 2 
An m X n feedback matrix, G, is defined as 

u = u c + Gx 

where u c is a command input. Then equation (34) becomes 


x = [F + BG]x + Bu c (35) 

The eigenvalues of the matrix F in equation (34) give the roots of the characteristic 
equation of the open-loop system, and the eigenvalues of the matrix [_F + BGJ in equa 
tion (35) give the roots of the characteristic equation of the closed-loop system. The 
CONTROL program uses the IBM Share program QREIG (i.e. , the QR algorithm) to 
fin d the eigenvalues. Setting 


G — aG-^ + bG£ 

where G and G 9 are m x n matrices and a,b - 0, 1,2 * the eigenvalues 

i 2 
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of [F + BG] yield root loci for the two feedback quantities associated with the matrices 
and G2- 
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TABLE 1. -APPROXIMATE OPEN-LOOP LITERAL TRANSFER FUNCTIONS 


1 K s K q / s 2 




( s) = -i- - KK s K q 


A(s) 12 M ( K + k 2 ) 


lKK s K q 


= -1- 

e A < s ) 12A 2 £ 2 (K + K 2 ) 


+ <^f S + 1 j( s2 + 2 ^t s + ^t 2 ) 


A(s) = (s +P 3 )(r s s + iJs—+ — £ s + 1^4 + +1 


, w c 


wg W2 


where 


2 K 2 
w t = T 


IK c 

£ ° ” 24a2;2 "c 



kk 2 k c 

Z 2 B v 

P 3 

12A 2 ^ 2 (k + K 2 ) 

21 

2 

co s = 

k + k 2 

I 

2 K 

W P ,2 

t rrij 

, ,2 - 

48/3 e A 2 ; 2 ^K + K 2 ) 

£ P 2w p\ 


IVK 

1 

“c 2 

1 1 
= — J + K 

w S 

2 VK h 
w 2 = 




B v / i _ z£\ 

2 ojp \nip I J 


\ 2 } 2 


IVi 2 m r 


^2 " 


VBv 

96/3 p A^ 
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TABLE 2. -RUDDER MODEL CONSTANTS 


K 

l 

V 

A 

K v 

K q 

K 1 

K 2 

K 


s 

l 6 r 

I 6 r 


1.8 slug-ft 2 
0. 166 slug 

10, 800 ft-lb/rad 

0. 158 ft 

* . 3 
5 m 

1.0 in 2 

50. 000 psi 

3 ~ 1 l ' 

3.05 in -sec /in. ^psi 

3 -1 

167 in -sec /in. 

940 ft-lb/rad 
1340 ft-lb/rad 
0. 003 in/V 
0. 0025 sec 

< 30° 

< 109 deg/ sec 


Nonlinear model 
B 500 lb/ft-sec -1 

V 

B c 40 lb 

K c 0. 0278 in 3 -sec Vpsi 


Linear model 
By 2000 lb/ft-sec * 

B c 0 

K c 0. 000084 in^-sec Vpsi 

x 0. 003 in. 

b max 
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TABLE 3. -COMPARISON OF APPROXIMATE LITERAL FACTORS AND EXACT FACTORS 


ADDroximate literal factors 


Case I 

Case n 

Case III 



Exact 

Approximate 

Exact 

Approximate 

Exact 

Approximate 

w s 

^K + k 2 

77.46 

77.46 

77.46 

77.46 

82. 12 

82. 12 

"h 

l 48 P e A 2 l 2 (K + K 9 ) 
yj IVK 




81.59 







86.5 


_i_ = _J__ _1_ 

ug u| + U] 2 

56. 18 

56. 18 

56.7 

56. 18 

63.0 

59.57 


K C I 

24 A 2 ; 2 i ° c 

.0139 

.0139 

.068 

.0139 

. 059 

. 0148 

P 3 

KK2K c 





.297 


12 A 2 ; 2 ( K + K 2 ) 





.331 

u?2 

|VK + 48/3 e A 2 Z 2 
V IW 2 m p 



2330 

1747 

2330 

1747 

h 

VBy ( 1 Vco 9 



2.58 




96^ e A 2 U + K 2 r ca, 2 



1.91 

2.58 

1.91 

“p 




1490 

1614 

1490 

1614 

“t 

IWj_ z 2 \ 

2w p \m p I / 

V?- 



4.02 

3.12 

4.02 

29.4 

3.12 

27.28 





h 

« 2 Bv rr 

21 > /k 2 


— 

— 

— 

.54 

.508 

z i 

K/3 

i 2 By 

"1 

*1 

K + K 2 



32.62 

21.14 

.0825 

.0775 
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Figure 1. Aircraft control -surface servo model. 



















P L = p r p 2 




Servo and load dynamics 



Figure 2. Closed-loop block diagram of an aircraft control -surface control system with position feedback 
and dynamic pressure feedback. Case I; m D = Ko = (3y = 0. 














Rudder surface 
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Figure 6. Rudder surface and electrohydraulic actuator. E -19830 
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Figure 7. Closeup of rudder electrohydraulic actuator. 











(c) — = 33.2 rad/ sec 
T p 


Figure 8. Concluded. 






Figure 9. Ground test frequency response of the rudder, the nonlinear model, and the linear model 
6 rc = 1.0°, peak to peak; K x = 1090 V/ft; K p = 0.0048 V/psi; = 16.6 rad/sec. 



Figure 9. Concluded. 




Figure 10. Ground test frequency response. 6r c = 1.0°, peak to peak; K x = 1090 V/ft; K p = 0.00785 V/psi 
— = 8. 3 rad/ sec. 




Figure 10. Concluded. 






(a) Simulated ground test. 


(b) Simulated flight test. 

Figure 12. Nonlinear model step response in simulated ground and flight tests. 
K x = 1090 v/ffc; Kn = 0. 0048 V/psi; — = 16.6 rad/see. 
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Figure 13. Continued. 





jw, rad/sec 'H rad/sec 



Figure 13. Concluded. 
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Figure 16. Concluded. 



Figure 17, Hinge-moment loading of rudder surface due to (3 response 
K v = 1090 V/ft; K d = 0.0048 V/psi; — = 16.6 rad/sec. 







Figure 18. Response of nonlinear model to simulated atmospheric turbulence. 

K x = 1090 V/ft; Kp = 0.0048 V/psi; — = 16.6 rad/sec; ov = 5 ft/sec; -rr£ = 10 sec. 
” Tp & V a 
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Figure 15, Concluded 






E 


Figure 16. Concluded. 


Figure 17. Hinge-moment loading of rudder surface due to f3 response 
K v = 1090 V/ft; K d = 0.0048 V/psi; — = 16.6 rad/sec. 
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Figure 18. Response of nonlinear model to simulated atmospheric turbulence. 

1 , L 

K v = 1090 V/ft; K n = 0.0048 V/psi; — = 16.6 rad/sec; a s = 5 ft/ sec; -rf- = 10 sec. 

A F Tr. & V a 
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